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ABSTRACT 

Observed chemical (anti)correlations in proton-capture elements among globular cluster stars are 
presently recognized as the signature of self-enrichment from now extinct, previous generations of 
stars. This defines the multiple population scenario. Since fluorine is also affected by proton captures, 
determining its abundance in globular clusters provides new and complementary clues regarding the 
nature of these previous generations, and supplies strong observational constraints to the chemical 
enrichment timescales. In this paper we present our results on near-infrared CRIRES spectroscopic 
observations of six cool giant stars in NGC 6656 (M22): the main objective is to derive the F content 
and its internal variation in this peculiar cluster, which exhibits significant changes in both light- and 
heavy- element abundances. We detected F variations across our sample beyond the measurement 
uncertainties and found that the F abundances are positively correlated with O and anticorrelated 
with Na, as expected according to the multiple population framework. Furthermore, our observations 
reveal an increase in the F content between the two different sub-groups, s-process rich and s-process 
poor, hosted within M22. The comparison with theoretical models suggests that asymptotic giant 
stars with masses between 4 and 5Mq are responsible for the observed chemical pattern, confirming 
evidence from previous works: the difference in age between the two sub-components in M22 must be 
not larger than a few hundreds Myr. 

Subject headings: globular clusters: individual (NGC 6656) — stars: abundances — :stars: AGB and 
post-AGB — stars: Population II — 



1. INTRODUCTION 

Although Galactic globular clusters (GCs) display a 
distribution in their global parameters (e.g., mass, metal- 
licity, concentration, horizontal-branch morphology) , the 
internal variation of elements affected by proton captures 
(here after p-capture elemen ts) appears a ubiquitous fea- 
ture (jCarretta et al.|[2C)09al b). It is clear that GC stars 
exhibit large changes in the C, N, O, Na, Mg, and Al 
abundances, whereas (in archetypical systems at least) 
internal spreads in iron-peak, heavy a- (Ca, Ti) and slow 
neutron-capture (s-process) elements all remain within 
observational uncertaint ies (Car retta et al . 2009c, 2010a; 
iJames et al.|[200l[S^^.200& : D'Orazi et al..,2010) . The 
changes in the p-capture elements give rise to a clear 
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chemical pattern: depletion in C, O, and Mg abundances 
always correspond to enhancements in N, Na, and Al (the 
so-called light-element anticorrelations) . This behavior 
bears the evidence of H burning at high temperature and 
points to the presence of multiple stellar generations. It 
is argued that the ejecta from a fraction of first gener- 
ation of stars (initially C-O-Mg rich, sharing the same 
chemical composition of field stars at the same metallic- 
ity) mix with primordial gas, providing a medium from 
which the second generation stars (C-O- Mg poor and N- 
Na-Al rich) formed. We refer the reader to lGratton et al.l 
l|20Tl for a recent, extensive review on this topic. In this 
scenario, the H-burning abundance patterns from the 
first generation stars are imprinted in the second genera- 
tion and are present from birth. The nature of the stars 
that enriched the intercluster gas remains uncertain but 
possible candidates include intermediate-mass asymp- 
totic giant branch (AGB) stars undergoing hot bottom 
burning (HBB, e.g..lD'Ant ona et al.l ll983: Ventura et alj 
(2001), fast rotating ma ssive stars (e. g., Dccrcssin et al] 
2007) , massive binari e s (Ide Mink et al . 2009), and novae 
((Smith fc Kraftlll996l: iMaccarone ~urek.2012i) . 

Interestingly, this already complex picture is fur- 
ther obfuscated by the pr esence of s ome peculiar clus- 
ters, such as u) Centauri (jJohnson fc Pilachowski 201(1 
Marino et all l2011aD NGC 1851 (lYong fc Grundahl 
2008t iCarretta et all [2010ai) . Terzan 5 (iFerraro et all 
2009h . NGC 671 5 (M54. ICarretta et al.l I2010bir 
and NGC 2419 (iCohen et all 120111 see however 
iMucciarelli et all 120121 for a different view). In these 
GCs, along with changes in p-capture elements, internal 
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variations in the heavy element abundances have been 
detected. Species ranging from the iron-peak (e.g., Fe) 
to the s-process elements (Ba, La) vary stochastically 
from cluster to cluster beyond what is expected from ob- 
servational errors. 

The metal -poor GC NGC 6656 (M22, [Fe/H]3~-1.70 
iHarrisI 119961 -updated in 2010) belongs to this class 
of GCs, and due to its peculiar nature h as re- 
ceived extensive attent io n (iPilachowski et"aL[ 
Norris fc F reeman 1983; 'Brown & Wallerstein' 
Kavser et al . 2008^ Marino et al. 2009: 
2009HDa Costa fc Mai±K)ll201C)i:lA!^es-Britoet al.ll201 



Recently, Marino et al. (2011b, hereafter MSKll) pre- 
sented results from their high-resolution spectroscopic 
study of 35 giant stars, deriving abundances for iron- 
peak, a, p-capture, and neutron-capture elements. This 
detailed abundance analysis provided a unique oppor- 
tunity to investigate the chemical enrichment history of 
M22. 

This GC is comprised of two distinct groups of stars, 
characterised by an offset in metallicity and in s-process 
element content. The first group displays a metallicity 
of <[Fe/H]>=-1.82±0.02 with <[s/Fe]>==-0.01±0.01 
and the second group has <[Fe/H]>=— 1.67±0.01 with 
<[s/Fe]>=-H0.35±0.02. Note that the [s/Fe] ratios were 
computed by averaging the abundances of Y, Zr, Ba, 
La, Nd (see MSKll for details). Each of these two 
sub-groups exhibits the classical Na-0 and C-N anticor- 
relations shown by the archetypical GCs. Given that 
the [Eu/Fe] ratio serves as a rapid neutron-capture (r- 
process) traceiQ and that enhancements in the s-process 
elements are not accompanied by a similar trend in the 
[Eu/Fe] ratio, it can be inferred that the dichotomy in 
the n-capture abundances may be due to a first genera- 
tion of polluters that produced the s-process only. One 
possible scenario, which was advocated by MSKll, is 
that the weak s-process component activated in stars 
with masses larger t han ^25 Mfr^ during core He-burning 
and C-shell phases (iRaiteri et al.l [19931 : iPignatari et al.l 
may have contributed to the observed abun- 
dance patterns . How ever, in a complimentary study, 
iRoederer et al.l (|2011[ ) focused on the heavy element con- 
tent (from Y to Th) of six stars across the two stellar 
sub-groups, and ruled out the massive star origin. They 
concluded that the gas from which the second stellar pop- 
ulation formed was enriched in s-process material from 
a class of relatively massive AGB stars (M w 5 Mq). In 
these stars, the production of s-process elements is due 
to the activation of the ^^Ne(a,n)^^Mg neutron source, 
whereas in their lower-mass counterparts t he main neu- 
tron source is the ^^C(a,n)^^0 reaction (jBusso et al.l 
[1999). 

Neither of the proposed scenarios provides a compre- 
hensive explanation for all the observed chemical features 
and we are left with numero us unsolved issues. For ex- 
ample, IRoederer et al.l ()2011f ) question why the s-process 
dichotomy is only present in M22: if massive AGB stars 
are the cause of the GC Na-0 anticorrelation, then all 
clusters should present the s-process elements correlated 

^ We use the notation [X/H]=A(X)-A(X)0, where 

A{X)=log{^)+12 

In t he Solar Syst em, Ri 97% of Eu was synthesized via the 
r-process I IBurris et al.. 200 0. and references therein). 



with Na an d anti-correlated wi th O, which is not ob- 
served (e.g.. lD'Orazi et al.|[20Tol) . 

In this paper we turn to an alternative diagnostic. We 
present fluorine abundances for a sample of six cool gi- 
ant stars in M22, carefully selected from both sub-stellar 
groups as defined by MSKll. Fluorine abundances are 
a powerful tracer of the polluter mass range in M22 be- 
cause the F production is highly dependent on the stellar 
mass. 

Th eoretical models of AGB stars (e.g., iJorissen et al.l 
119921 ) predict that F is produced due to the activation 
of the chain of reactions "'^^0(p,q;)^^N(q;,7)^^F in the He 
intershell during the recurrent thermal pulses associated 
with He burning. During the early phases of each ther- 
mal pulse, the H-burning ashes are ingested in the con- 
vective region developing in the He intershell. These 
ashes are rich in ^'^C and ^^N and their ingestion in the 
He- burning layer results in production of due to a 
captures on ^^N. At the same time protons are released 
by the ^^N(n,p)^*C reaction, with neutrons coming from 
^'^C(a,n)^^0. After the quenching of each thermal pulse 
the envelope may sink in mass deep in the He intershell 
and carry ^^F to the convective envelope via a process 
known as the "third dredge-up" (TDU). The peak of F 
production in AG B stars is reached for stars of initial 
masses ~ 2 Mq (|Lugaro et al.l I2004D . If the mass of 
the stars is higher than roughly 5 Mq, and depending 
on the metallicity, fiuorine is destroyed both via a cap- 
tures in the He intershell, and via proton captures at the 
base of the convective envelope due to HBB. AGB stars 
that experience HBB can also destroy O and Mg and 
produce Na and Al; as a consequence, according to the 
multiple population scenario, we should expect the abun- 
dances of F to be correlated with O (and Mg) and anti- 
correlated with those of Na (and Al). This predict ion 
was observationally confirmed by iSmith et al.l (I2005D in 
the in termediate-metallicity GC M4 and bv lYong et al.l 
(|2008f) in NGC 6712. 

Other sites have also been suggested for F pro- 
duction: the ^^-pr ocess in core-collapse supernovae 
(iWooslev et al.lll990D. and core H e burning in Wol f -Raye t 
stars (|Mevnet fc Arnouldl 120001: iPalacios et all 120051 ). 
However, low-mass AG B stars are the only site ob- 
serva tionally confirmed porissen et al.|[T993 lAbia et all 

[Moh . 

This paper is organized as follows: observations are de- 
scribed in Section [21 while details on abundance analyses 
are given in Section l3l Our results are then presented in 
Section [J] and discussed in Section [5l A summary closes 
the manuscript (Section [H]). 

2. OBSERVATIONS 

Our sample includes six giant stars, for which stellar 
parameters (Tos, logg, [Fe/H] and microturbulence ^) 
along with p-capture and s-process element abundances 
were derived by MSKll. We selected three stars be- 
longing to the metal-poor component (MP, also s-process 
poor) and three stars from the metal-rich (MR, s-process 
rich) one. Within each of these sub-groups we also se- 
lected both 0-rich (Na-poor) and 0-poor (Na-rich) stars, 
spanning a range from [O/Fe]=-t-0.11 to [0/Fe] =-1-0.48 
dex. 

The main objective of our investigation was the de- 
termination of the fiuorine abundances. Observationally, 
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Fig. 1. — Example of a spoctmm for star III-52 

F (whose only stable isotope is ^^F) is difficult to detect 
spectroscopically; the only atomic lines (the ground-state 
transitions of F i) that might be revealed lie in the far 
UV. On the other hand, HF molecular transitions are 
easily observable in the near-infrared (around ~23000 
A). Our analysis focuses on the HF(1— 0) R9 transition, 
located at A23358.3 A. Although not being the strongest 
feature in this wavelength range, this line is considered to 
be one of the be st abundance ind i cator for F, because i t 
is free of blends ()Abia et al.ll2009l: iLucatello et"aIll201lD . 

High-resolution, near-infrared spectroscopic observa- 
tions were carried out in service mode with CRIRES 
(CRyo genic high-re solution InfraRed Echelle Spectro- 
graph, iKaeufletaDHOQl) located at VLT UTl on 2011 
April, July and August (program: 087.0319(A), PI: VD). 
In Table [U we list information on target stars, report- 
ing identifications, magnitudes (see MSKll for details), 
exposure times, and S/N ratios per pixel around the HF 
feature. Note that our selection was limited to the cooler 
stars in the MSKll's sample due to our imposed require- 
ment of relatively strong HF lines. We also observed 
several (hot) early-type stars before and/or after of each 
target observations, in order to remove telluric contami- 
nation from our scientific frames. 

We employed the 0.4 slit and the grating order #24, 
achieving a resolution of i?~50,000 and a wavelength cov- 
erage from A22948.5 A to A23410.3 A. This allowed us 
to include the HF(1— 0) R9 line, numerous ^^C^^O vibra- 
tion/rotation lines (used to derive C abundances) and the 
Na I line at A23379 A. 

Data reduction was accomplished by means of 
the CRIRES pipeline (version 2.4), running under 
the GASGANO environmenlF^. which provides one- 
dimensional, wavelength calibrated spectra. Telluric fea- 
ture subtraction, rest-frame translation and continuum 
normalization were then carried out within IRAlQ. An 

http:/ /www. eso.org/sci/software/gasgano/ 
IRAF is the Image Reduction and Analysis Facility, a general 
purpose software system for the reduction and analysis of astro- 



example of our spectra is shown in Figure [T] for star III- 
52; HF, CO and Na i lines are marked. 



3. ABUNDANCE ANALYSIS 

Fluorine abundances were determined through spec- 
tral synthesis using the MOOG code (Sncden 1973, 2011 
version) and the Kurucz (1993) set of stellar atmosphere 
models (with no overshooting) as in MSKll's analysis, 
from which we retrieved stellar parameters as well as O 
and s-process element abundances. How ever, had we in- 
stead adopted the MARCS grid (Gustafs son et al.| [2008). 
the difference on the resulting abundances would have 
been less than 0.04 dex. 

Concerning the HF feature, we took as excita- 
tion potential (EP) y=0.227 eV (iDecinl I2000D and a 
log5/=-3.971 (jLucatello et al.ll20lH ). This last value is 
very close to that used in previous stu dies focu sing on F 
abundance deter mination in GCs (e.g.,lSmithet al. 2005]; 
lYong et al.' '2008: Alves- Brito etaLll2012r ). logqf =-3.955 
(from Ijorissen et al.,.199^ . On the other hand, our EP 
is about 0.25 eV lower compared to those works: this 
difference implies an offset between our abundances and 
those previously published in the literature of roughly 
0.30 dex (see also the discussion in Section [4.11) . Fur- 
thermore, we derived C abundances. For this purpose, 
we assumed O values from the optical range given by 
MSKll, since our spectra did not cover any suitable OH 
line, whose stronger transitions extend in the H band 
(around -15000A). The CO line lists come from B. Plez 
(private communication). Finally, for the Na i line at 
A23379, atomic parameters (x=3.750 eV; log5/=0.530) 
were taken from VALEH. 

As a first step, we checked our lin e list on the infra- 
red atlas of the Arcturus spectrum (Hinkle et al."1995|, 
available at ftp:/ /ftp. noao.edu/catalogs/arcturusatlas/). 
Assuming a Teff=4286 K, log.(7=1.67, ^=1.74 kms^^ 
and [Fe/H]=— 0.52 (following iRamirez fc Allende Prietol 
1201 If ) we obtained an [F/Fe]=— . 1^3, to be compared 
to the value given bv lAbia et all (|2C)09( ) of [F/Fe]=0.10 
dex. The difference is completely explained by the higher 
EP adopted in that study. 

Moreover, we inferred a C abundance of A(C)=8.01 
(under the assumption that A(0)=8.81, that is 
[O/Fe]=0.40, with a solar abundance of A(0) 0=8.93), 
which is i n ver y good agreement with values of 
lAbia et al.l (pOOl (i.e., A(C)=8.06) and IRvde et al.l 
(|2010D (A(C)=8.08). As for C solar abundance 
we adopted the value of A(C)0=8.56, leading to 
[C/Fe]=-0.03. 

Finally, from the Na i feature at 23379A, we derived 
A(Na)=6.01 which resuhs in [Na/Fe]=+0.2 dex (setting 
A(Na)Q=6.33). 

Comparison between synthetic and observed spectra 
were carried out in a similar way for our sample stars; 
an example of spectral synthesis is given in Figures [2] 
and [3] for star C. 

nomical data. IRAF is written and supported by National Optical 
Astronomy Observatories (NOAO) in Tucson, Arizona. 

Vienna Atomic Line Database 

(www.astro.uu.se/ vald/php/vald.php) 

^"^ We used as solar fluorine abunda nces A(F)0=4.56 
HAnders fc Grevess6|[T989l : lAspIund et~ani200a) . 
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TABLE 1 

Information on Target Stars 



Star 


RA 


Dec 


V 


K 


Exposures 


S/N 




(hh:mm;ss) 


{"■'■■") 


(mag) 


(mag) 


(s) 




IV-97 


18:36:41.06 


-23:58:18.9 


11.043 


6.759 


10x90s 


300 


III-14 


18:36:15.10 


-23:54:54.6 


11.134 


6.743 


17xl20s 


400 


III-15 


18:36:15.61 


-23:55:01.2 


11.362 


7.138 


8xl80s 


300 


C 


18:36:10.21 


-23:48:44.0 


11.309 


6.737 


6xl80s 


300 


III-52 


18:36:10.18 


-23:54:21.8 


11.526 


7.459 


10x180s 


350 


V-2 


18:36:28.02 


-23:55:01.6 


11.498 


7.276 


14xl20s 


350 




A(F) = 1.92; 2.07; 2.22 





1 


1 


1 






y 

1 




\ 














A(C) 


= 6.07, 6.27, 6.47 

1 .... 1 


1 





23358 

Wavelength (A) 



23300 23310 23320 23330 

Wavelength (A) 



Fig. 2. — Synthesis of the HF feature for star C. 



The sensitivity of the F abundance to input stehar pa- 
rameters was evaluated by separately changing effective 
temperature, surface gravity and microturbulence values. 
The intensity of the synthetic line of the HF is particu- 
larly sensitive to the adopted Tcft, the other parameters 
affecting it at a lower degree (see also lAbia et all 120091 : 
2011). A change of A(Tcff )=-l-70 K, A(log5)=-h0.15, and 
A(^)=-|-0.13 km s~^ (conforming to error estimates given 
in Table 4 of MSKll) results in a difference in A(F) of 
+0.10, 0.02 and —0.02 dex, respectively. The variation 
of the input metallicity in the model atmosphere has in- 
stead a negligible effect. These are the typical uncertain- 
ties that we then summed in quadrature providing a total 
error, due to stellar parameters, of 0.11 dex in our [F/H] 
ratios. Errors due to the best-fit determination (related 
to the S/N of the spectra and including uncertainties 
due to the continuum placement) are instead ±0.07 dex. 
However, we caution the reader that this value should be 
treated as a lower limit, since the impact of the telluric 
correction on this considerably weak feature is significant 
(see Section |4|. 




23376 2337B 23360 23362 

Wavelength (A) 

Fig. 3. — Same as Figure [2] but for CO (upper panel) and Na 
(lower panel) 
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Fig. 4. — Fluorine abundances ([F/H]) as a function of [O/H]. 

The total internal error in [F/H] is then obtained 
adding in quadrature both uncertainties, resulting in 0.13 
dex. 

Finally, as far as C and Na are concerned, the typical 
micertainties are 0.12 and 0.10 dex, respectively. 

4. RESULTS 

Our results are shown in Table [21 where we report stel- 
lar parameters and abundances from MSKll along with 
our estimates for F, C, and Na. Even within our quite 
limited sample (six stars), we found that the F abun- 
dance shows a large star-to-star variation, ranging from 
[F/H]=-2.82 dex to [F/H]=-2.23 dex (i.e., a factor of 
~4). The average abundance is <[F/H]>=-2.55±0.08 
(rms=0.20), implying that the amplitude of this change 
is beyond the measurement uncertainties (see Section [3]). 
Moreover, taking into account the typical errors for O 
and F abundances, our study suggests that the F vari- 
ation is compara ble with that of O, as also found by 
lYong et all ((200l) in the GC NGC 6712. 



In Figures m and [S] we plot [F/H] ratios as a function of 
[O/H] and [Na/H], respectively. As one can see, F abun- 
dances are positively correlated with O: considering the 
whole sample, the Pearson correlation coefficient results 
in r=0.89, with a probability to be random smaller than 
2%. Focusing on the F-Na diagram (Figure [5]), there 
is the hint for a F-Na anticorrelation, but in our small 
sample the linear correlation coefficient is not statisti- 
cally meaningful. However, this does not prove the lack 
of an anticorrelation, because we are dealing with small 
numbers (only six points) , heavily reducing the power of 
statistical tests. In addition, there is no a priori reason 
why we should combine the two sub-groups as far as the 
F-Na plane is concerned, because we do not expect that 
they must behave in the same way. If we look at each 
component separately the presence of a F-Na anticorre- 
lation is indeed much more evident (given that we have 




-2.4 -2 -1.6 -1.2 -O.E 



Fig. 5.— [F/H] versus [Na/H] from MSKll (left panel) and from 
this study (right panel). Symbols are as in Figure |4] 

only three points for each group, it is not meaningful to 
perform a statistical test separately). 

The observed chemical pattern can be explained as the 
evidence of H-burning at high temperature, via the CNO 
cycle, which causes the destruction of F, in conjunction 
with O depletion and Na enhancement. IVIore interest- 
ingly, those (anti)correlations are revealed in each of the 
]VI22 sub-components (the s-rich and s-poor groups); the 
implications of this finding are discussed in detail in Sec- 
tion [5l 

Regarding Na, we show both our estimate from near- 
infrared spectroscopy as well as LTE abundances from 
the optical range by IVlSKll (right and left-hand panels 
of Figure E]). A difference of A([Na/Fe])=0.31±0.04 
(rnis=0.09) dex is found between the two estimates 
(see Figure [6]) where we compare the two mea- 
surements; non-LTE effe cts can totally account for 
such a discrepancy (e.g., iLind et all 120 111) . The to- 
tal average Na abundance is <[Na/Fe]>=0.05±0.13 
(rms=0.33); considering separately the two groups 
we obtain instead a <[Na/Fe]>s_poor=— 0.12±0.23 
and <[Na/Fe]>s-rich=+0.22±0.08, which implies 
A™;^,[Na/Fe]=0.34±0.17 dex. This value is in 
good agreement with that derived by JVISKll, 
based on the whole sample of 35 giants, being 
^poor[Na/Fe]=0.23±0.07. 



Finally, our sample displays an average 
C abundance of <[C/Fe]>=-0.74±0.18, 
with <[C/Fe]>^_poor=-l-12±0.02 and 

<[C/Fe]>s_rich=-0.37±0.14. Thus, on average, 
s-rich stars exhibit larger C abundances, which quali- 
tativ ely agrees with previous works (e.g., IBrown et al.l 
119901 ). However, while the difference between the two 
groups from the optical CH bands derived by JVISKll 
is A™J;,[C/Fe]=0.35±0.13 dex, we achieved a much 

larger value of Agg^JC/Fe]= 75ib0.10 dex. Th e same 
conclusion was reached by lAlves-Brito et al.l (|2012j ) 
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TABLE 2 

Stellar Parameters and Elemental Abundances 



Star 


Toff* 


log 3* 


[Fe/H]* 


€* 


[O/Fe]* 


s-rich* 


[C/Fe] 


[F/Fe] 


[Na/Fc] 




(K) 


(cms^^) 




(kms^-'-) 












IV-97 


4000 


0.05 


-1.94 


2.00 


0.40 


no 


-1.10 


-0.70 


-0.40 


III-14 


4030 


0.35 


-1.82 


2.15 


0.48 


no 


-1.10 


-0.80 


-0.30 


III-15 


4070 


0.40 


-1.82 


1.85 


0.11 


no 


-1.15 


-1.00 


0.34 


C 


3960 


0.30 


-1.69 


2.25 


0.25 


yes 


-0.60 


-0.80 


0.30 


III-52 


4075 


0.60 


-1.63 


1.75 


0.45 


yes 


-0.10 


-0.60 


0.05 


V-2 


4130 


0.65 


-1.57 


1.75 


0.15 


yes 


-0.40 


-0.90 


0.31 



* from MSKll. 




-0.5 5 



[Na/Fe]„,,,, 

Fig. 6. — [Na/Fe] from the optical range by MSKll and from 
this study 

who found a variation of A™|;jC/Fe]=0.78±0.15 dex, 
from high-resolution, near-infrared spectroscopy of nine 
cool giants (see also Section I4.ip . The reason of such 
a difference in C abundance from the optical and from 
the near-infrared is not clear and no obvious trends 
with stellar parameters (e.g., temperature, gravity, 
microturbulence, and/or metallicity) seem to be present. 
Further investigations are needed to explore this issue. 

4.1. Comparison with \Alves- Brito et al\ !l201i ) 

In a recent paper, lAlves-Brito et all (j2012D carried 
out high-resolution (i?=50,000), near-infrared (both 
H and K bands) spectroscopic observations with 
Phoemx@Gemini-South of nine red giant branch (RGB) 
stars in M22. They investigated their F content, present- 
ing also C, N, O, Na and Fe abundances. 

Four of our six stars are in common with that study, 
namely, III-14, III-15, III-52 and IV-97. For stars 
III-14 and III-52, those authors inferred [F/Fe]=-0.40 
and [F/Fe]=— 0.20 dex, respectively, while we obtained 
[F/Fe]=-0.80 and [F/Fe]=-0.60 dex. The adopted EP 
value can easily justify this divergence, accounting for 
about 0.30 dex (see Section[3l); the source of the remain- 



ing '^0.1 dex can be ascribed to the continuum place- 
ment, which is critical in determining abundances from 
such a w eak line. On the other hand, for stars III-15 
and IV-97 lAlves-Brito et al.l(|2012{ ) obtained [F/Fc]=0.28 
and [F/Fe]=0.25 dex, entailing discrepancies with our 
estimates larger than a factor 10. Note that stellar pa- 
rameters (Toff, logg, and ^) are the same in both works, 
as they come from the analysis of MSKll; for the in- 
put metallicity, Alves-Brito et al. used instead their own 
values coming from the IR spectroscopy and showing an 
offset of -1-0.13 dex compared to the optical ones. We 
investigated the nature of this substantial discordance, 
and attributed it to the telluric feature subtraction. In 
the upper panel of Figure [71 we directly compare our 
spectrum for star III-15 (solid line) with that used by 
Alves-Brito et al. (dotted line). Their spectrum presents 
stronger features in the vicinity of HF line, features ex- 
pected from telluric contribution. To completely remove 
the contamination, Alves-Brito et al. realized that they 
needed early-type star targets before and/or after each 
scientific frame but the logistics of their run made this 
very difficult. We instead could observe such targets: 
the correction to our data, in turn affects the placing 
of the continuum and removes many strong features (as 
shown in the lower panel of Figure [T]). As expected, 
such an effect is significant for the HF feature, due to 
its intrinsic weakness, but only marginally affects the C 
and Na abundance determinations (due to the strength 
of their lines), and hence most of the conclusions of that 
paper. This is shown in Figure El where we plot our 
[X/Fe ] ratio as a function of those from lAlves-Brito et al.l 
()2012f ) for the four stars in common: C and Na are 
comparable between the two studies, with differences 
of A[C/Fe]=-H0.20±0.17dexand A[Na/Fe]=-h0.18±0.11 
dex (in the sense Alves-Brito's study minus our val- 
ues); if we take into account the offsets in [Fe/H] 
and in the adopted solar abundances (they assumed 
A(C)0=8.42 and A(Na)0=6.17), those values become 
A[C/Fe]=-f'0.19±0.17 and A[Na/Fe]=-H0.15±0.11 dex. 
On the other hand, discrepancies in F are signifi- 
cant and cannot be recovered from the different EP 
and/or solar abundances, being for the whole sample of 
A[F/Fe]=-|-0.75±0.19 dex. 



Fluorine variations in M22 



7 



T 




,8 I . , , , , , . . , , , , , . 1 , , , , 1 

23345 23350 23355 23360 23365 



Wavelength (A) 

Fig. 7. — Upper panel: comparison of our spectrum for star III-15 
(solid line) with the one by Alvcs-Brito ct al. (201^) (dotted line). 
Lower panel: superimposition of our spectra for the same sample 
star with and without the telluric line subtraction (dot-dashed and 
solid lines, respectively). 




Fig. 8. — Comparison for the four stars in common with Alves- 
Brito et al. (2012, here labeled as AB12). Symbols are as in 
FiguresHEllg 

5. DISCUSSION 

Our main result is the detection of F variations across 
our sample significantly larger than the observational un- 
certainties. As shown in Figures |3] and [5l the changes 
in the F abundances are correlated with O and anti- 
correlated with Na. This chemical pattern qualitatively 
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matcnes tne predictions from tne multiple population 
scenario, according to which the stellar ejecta from which 
second generation stars formed carry the signature of hot 
H burning causing enhancements in Na (N and Al) and 
depletions in O and F (C and Mg). The F-0 diagram 
presented in Figure [S] demonstrates that M22 shares a 
similar behaviour as M4 and NGC6712, the other two 
GCs for which F has been explorecpl. 

Furthermore, the F-O-Na (anti) correlations can be 
marked separately within each of the two sub-groups en- 
closed in M22 as clearly illustrated in Figures |4] and \5\ 
where the s-process poor and the s-process rich stars are 
labeled with empty and filled symbols, respectively. The 
same conclusion was drawn by MSKll when considering 
the Na-0 and C-N planes. 



Very interestingly, beyond the internal spread in F 
characterizing each sub-component, we measured an in- 
crease in the F content between the two different stel- 
lar generations in M22. The s -process-rich group has, 
on average, larger F abundances than the s -process-poor 
group. This is shown in Figure I10[ where we plot our F 
abundances ([F/H]) as a function of [La/H] from MSKll. 
There is a positive correlation between the two ratios, 
suggesting that the polluters responsible for the s-process 
production must account for a simultaneous F produc- 
tion. 

There are two classes of objects producing both s- 
process elements and fluorine. The first one is very mas- 
sive stars (mass roughly > 4OM0). These produce F 

ICunha et al.l II2003I ) presented F abundances for two giants 
in ui Cen. However, they provide an F measurement only for star 
ROA219, giving an upper limit for star ROA 324. Discussion re- 
lated to the internal F variation in this peculiar GC is still not 
possible with the currently available measurements. For this rea- 
son we acquired CRIRES spectra of 12 ui Cen giants; results will 
be presented in a forthcoming paper (S. Lucatello et al., in prepa- 
ration). 
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Fig. 10.— [F/H] versus [La/H]. Symbols as in Figures g] [5] H 
and|8] The dashed Hue is a least-squares fit to data points. 

in the initial phases of core He burning and expel it in 
the int erstellar medium via wind s during the Wolf-Rayet 
phase (jMevnet fc ArnouTdll2000[ ). They also produce s- 
process elements during core He and shell C burning 
(e.g., iPignatari et al. 2010). Production of both F and 
s-process elements in these massive stars depends on the 
initial CNO abundances and thus decreases with the stel- 
lar metallicity. Inclusion of stellar rotati on enhances the 
s-process production at low metallicity (jPignata ri et al.l 
[2OO8; Chiappini ct al. 2011); howe ver, it appears to de - 
crease the production of fluorine (iPalacios et al.l I2005D . 
One problem already stressed bv lRoederer et all ()2011[ ) 
when considering these stars is that there is no rea- 
son why the SNe that enriched the s-process rich group 
host the weak component and those that polluted the 
s-process poor do not (see Roederer et al. for details). 

The second class of objects produci ng both F and 
s-pro ces s elements a re AG B stars ([Forestini et al.l 
1992t IJorissen et al.l 119921: iMowlavi et al I 119981: 
Karakas fc Lattanziol 120031: iLugaro et al.l |2004 2012). 
To get deeper insights into the nature of the candi- 
date AGB stars possibly responsible for the observed 
abundance trends in M22 we compare our res ults with 
the recent set of models bv ILugaro et al.l ()2012.) . They 
presented AGB models for masses 0.9—6.0 Mq and 
metallicity three times lower than that of the cluster 
under consideration (i.e., [Fe/H]=— 2.3 dex). 

From our data we infer that there is an increase of 
A[F/H]^:i;;j,^„';=-h0.40±0.15 dex in the fluorine content 
between the two groups. This estimate was done by 
taking into account the F content of the 0-rich stars 
only, because they do not show any depletion due to 
the HBB. We averaged the F abundances in stars IV- 
97 and III-14, both belonging to the s-poor group, find- 
ing a mean value <[F/H]>=— 2.63±0.01 dex. Since we 
have only one 0-rich star in the s-rich group, IIT52, 
we chose it as representative of the F abundance for 
the group, that is [F/H]=-2.23±0.15 dex. The F in- 
crease of -1-0.40 dex is accompanied by a correspond- 



ing enhancement in La of A[La/H]^_™'j.=0.56 ±0.18 
dex, since the values are <[La/H]>=— 1.92±0.15 dex 
and [La/H]=— 1.36±0.10 dex, respectively, for the two 
groups. 

Comparing th e se va lues with the model predictions 
bv ILugaro et al.l (j2012l ) we found that AGB stars with 
masses of ~ 4—5 can well reproduce the observed 
pattern. Lower-mass AGB models do not fit our observa- 
tional requirements because they over-produce fiuorine. 
This is true even if we consider that these AGB model 
predictions are roughly one dex too high to match the ob- 
se rvation of carbon-enha nced metal-poor (CEMP) stars 
bv lLucatello et all ()2011| ). Production of F in the 4-5 
Mq mass range depends on the delicate balance between 
the operation of the TDU and of HBB. These stars suf- 
fer HBB and destroy F during the early phases of their 
AGB evolution, however, toward the end of the evolution, 
as the mass of the convective envelope decreases, HBB 
ceases while the TDU is still active resulting in mild F 
enhancements at the stellar surface. This explains why 
while the most prolific AGB stars in terms of fluorine 
production have initial masses around 2 Mq, F produc- 
tion still occurs at slightly higher masses. On the other 
hand, more massive AGB models (>5 Mq) experience 
hotter HBB and thus more efficient F destruction as well 
as higher temperatures in the thermal pulses activating 
also ^^F(Q;,p)^^Ne reactions. This, combined with fewer 
final TDU episodes when HBB has ceased means that 
they do not replenish F at the stellar surface. 

Interestingly, the same conclusion is drawn by 
IRoederer et al.l (j2011l ) by exploring the heavy-element ra- 
tios, [hs/ls] and [Pb/Zisj. Comparing their abundances 
with models by Roederer et al. (2010), these authors de- 
duced that the low mass AGBs (<3M0) cannot account 
for the observed trend. More importantly, they con- 
cluded that a match to the s-process element abundances 
is provided by the 5 Mq AGB m odel, and this iscon- 
firmcd by checking the models of ILugaro et al.l (|2012D . 
Furthermore, and very interestingly, these models pre- 
dict Na and C production in agreement with the obser- 
vations. 

Our result provides a further, independent confirma- 
tion to this previous hint: indications from both light 
(fiuorine here) and heavy elements converge toward the 
AGB stars of the same mass as the best candidate pol- 
luters, indicating that the age difference between the two 
sub-groups in M22 cannot be larger than a few hundreds 
Myr. It should be mentioned that by anal yzing the dou- 
ble su b-giant branch (SGB) of this cluster, IMarino et al.l 
(|2012l ) concluded that the age spread can be at most 
-300 Myr. 

The fact that three independent studies, involving dif- 
ferent and complementary techniques/approaches, pro- 
duce the same result is encouraging. However, a compre- 
hensive understanding of the whole picture is s till miss- 
ing. As also stressed bv IRoederer et all ()201lD . if rela- 
tively massive AGBs (4-5 Mq) produced the s-process 
elements, and if these stars are also responsible for the 
observed p-capture element anticorrelations, then it is 
not explained that an s-process enrichment is present in 
M22 but is not associated with O and Na abundance 
anomalies, nor it is seen in any other GCs where the 
Na-0 anticorrelations are observed. 
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On the other hand, it is clear from the observed F 
abundance trends that we must select slightly more mas- 
sive AGBs, i.e., > 6M0 if we wish to explain the light 
element variations in GCs, since this is the AGB mass 
range where F and O can be destroyed by HBB resulting 
in the observed F-0 correlation. We might tentatively 
suggest that perhaps the production of the s-process el- 
ements in AGB models with initial mass > 6Mq is less 
efficient than currently predicted. This could be the re- 
sult of a stronger mass-loss rate or a less efficient TDU 
in this mass range tha n those employed in the models 
by iLugaro et al.l (j2012( ) . This possibility is well within 
model uncertainties and need to be investigated; in a 
forthcoming paper (D'Orazi et al., in preparation) we 
will attack these issues, presenting new observations and 
AGB models and discussing their strength/ weakness in 
reproducing the observed abundance trends in GCs. 

Alternatively, we may conclude that massive AGBs 
are not the inter-cluster polluters; however several lines 
of evidence point to those stars, such as the need for 
a Li production between firs t and second generation 
stars (as in the case of M4, 'D'Orazi & Marino| [20T0I: 
[Muccia relli et al. 2011; Monaco et al. 20l"|). Further ef- 
forts, from both observational and theoretical perspec- 
tives, are needed; in particular the lack of a complete 
set of models for AGB stars and massive rotating stars 
with different mass and metallicity, following the whole 
nucleosynthetic path from Li to Pb, still hampers a quan- 
titatively robust comparison between theory and obser- 
vations. 

6. SUMMARY AND CONCLUDING REMARKS 
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